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ABSTRACT 


In this report, we discuss the physics of fading and 
summarize the methods to mitigate its devastating effects for 
communication channels. We evaluated performance curves of 
Differentially Encoded Binary Phase Shift Keying (DBPSK), Rate 
1/2 Convolutionally Encoded Differential Quadrature Phase 
Shift Keying (DQPSK) and Minimum Shift Keying (MSK), by 
Simulation for both AWGN and fading channels. Further, we 
Simulated typical transponder components of a UHF satellite 
System and evaluated their effect for the above modulation 
structures. Finally, we calculated the carrier-to-noise ratios 
by using simulation curves to assess link margins for the UHF 


Follow-On Satellite. 
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I. INTRODUCTION 


During the period when the HF, or so-called short-wave 
radio was the corner stone of long distance communications, 
Signal fading, also termed signal scintillations, was accepted 
aS an inevitable problem that one should live with. 

This problem was thought to be solved when the Early Bird 
(INTELSAT I) was launched in April 1965 by International 
Telecommunications Satellite Organization (INTELSAT) [Ref. 1: 
eee). 

Satellites, with their broad coverage and line-of-sight 
communications capabilities, enabled the communications 
community to utilize higher carrier frequencies. Ionospheric 
disturbances were thought to be negligible at these higher 
frequencies and were modeled using only additive white 
Gaussian noise (AWGN). But deployment of these earth-orbiting 
devices has generated reports of unexpected disturbances in 
the communication channels [Ref. 2: p. 5]. Because of those 
reports, some measurements were made all over the world, 
revealing that the cause for these disturbances was, in fact, 
1onospheric irregularities. This knowledge brought a new 
understanding and awareness for the importance of fading for 


UHF satellite communication systems. Today the importance of 


fading is well appreciated by the members of the 


COMMUNICaAETONS Conmunsey- 


A. FADING AND ITS EFFECTS 

Before getting into the effects of fading for satellite 
channels and the methods to mitigate these effects, we will 
try to summarize the physics and causative components of 
fading. 

Signal fading, or what radio astronomers call 
scintillation, is caused by signal propagation through a 
region which has a medium with random fluctuating index of 
refraction [Ref. 2: pp. 6]. These fluctuations result from 
atmospheric turbulance, meteorological events, dust along 
tropospheric paths and the ionization structure along 
transionospheric paths. These disturbances have two main 
effects over an electromagnetic signal, one of which is a 
random phase disturbance due to variations in the phase 
velocity and the other is a random magnitude fluctuation due 
to scattering. 

Recovering the transmitted information from this disturbed 
Signal is the main challenge for communications engineers who 
are in search of a fast and reliable communication system. If 
a communications system is not designed to operate in the 
presence of signal fading, the performance of the system may 
be severely degraded if fading is encountered. Some satellite 


communications systems were originally designed for operation 


in the AWGN channel based on the assumption that the 
ionosphere would have negligible effect on carriers above 100 
Mhz or so, but later UHF high frequency carrier systems were 
observed to be significantly affected by signal fading. 

Fading can be observed in two different forms related to 
the symbol period. Slow fading, in which the symbol amplitude 
and phase variations can be observed to slowly change compared 
to the symbol period, has a highly random fluctuating error 
pattern with long bursts of errors scattered among prolonged 
error free periods. Fast fading, with its rapid fluctuations 
during the symbol period, gives slight chance for successful 
coherent demodulation and the bit error rate approaches 0.5 
feet. 2: p. 7). 

Since the motion of the ions and index of refraction 
Changes are random to the channel user, these random 
variations can be treated in statistical terms. Fading 
Statistics also can be encountered in two different forms 
according to the physical environment of the propagation 
medium. A Rayleigh fading channel model is the most common 
multipath propagation model. In this model, signal fading is 
caused by time-variant multipath characteristics of the 
channel, thus forming a channel impulse response which can be 
modeled as a zero mean complex-valued Gaussian process with a 
Rayleigh distributed envelope [Ref. 2: p. 7]. When the medium 
has fixed scatterers or signal reflectors in addition to the 


randomly moving scatterers, the channel can no longer be 


modeled as having a zero mean. The resultant channel is the 
Ricean distributed fading channel. 

Also important is frequency selectivity. If the coherence 
bandwidth of the channel is very small in comparison to the 
bandwidth of the transmitted signal, the channel is said to be 
frequency selective [Ref.3: pp. 707-708). Otherwise, it is 
non-frequency selective. When fading is observed in the 
frequency selective form, the distortion of the signal is very 


severe compared to the nonselective form. 


B. FADING MITIGATION METHODS 

There are two major methods for mitigating the damaging 
and undesirable effects of fading on communications systems. 

1. Modulation Scheme Selection 

The first solution is directly related to the modulation/ 
demodulation method. With selection of a modulation scheme 
such as differentially coherent phase shift keying (DPSK), 
which does not require carrier phase tracking, instead of 
conventional coherent phase shift keying (CPSK), which is an 
optimum modulation scheme for AWGN channel, communications 
systems are less suspectable to the rapid phase disturbances. 
When a modulation scheme is selected in which the channel bit 
rate does not exceed the minimum value of the frequency 
selective bandwidth, the system is said to be working in its 
Slow fading limit [Ref. 2: p. 17 ]. This arrangement in the 


modulation/demodulation scheme improves performance to a limit 


which still may not be acceptable with normal link power 
Margins of the system, but is better than the fast fading 
performance of CPSK. In this event, a second mitigation step 
must be invoked. 
2. Use of Error Correction Coding and Interleaver 

This second step incorporates error correction coding 
and diversity techniques to improve average error rates ina 
slow fading environment. 

a. Error Correcting Coding 

The basics of error correcting coding were first 

introduced by Shannon in 1948. Since then, a great deal of 
effort has been expended on the problem of devising efficient 
encoding and decoding methods for error control due to channel 
noise. 

The objective of coding is to utilize redundant 
data to retrieve information during the decoding phase without 
increasing the average power. This is done by increasing the 
effective minimum distance between the Signals. One way to do 
this is to use convolutional codes. Convolutional codes have 
proved useful in mitigating effects of fading particularly 
when combined with interleaving or other forms of diversity 
techniques. When provided with a reasonable error pattern, 
these codes, in conjunction with soft decision Viterbi 
decoding, are effective in reducing binary error rates from 


10* to 10° or less. In Chapter II we provide the detailed 


specification of the rate 1/2 convolutional encoder we 
implemented. 
b. Diversity Techniques 

Diversity techniques use the concept of 
redistributing extensively long error bursts to present a more 
random-appearing error pattern to the decoder. These 
techniques split modulator output in a controlled manner and 
recombine after the channel to break the effect of burst 
fading patterns into shorter groups of errors for the decoder 
to accommodate [Ref. 2: p. 208 ]. Diversity techniques are 
effective for the slow fading environment in which the fading 
pattern can be observed in long burst forms. 

The block interleaving method useS aé_e two- 
dimensional memory array. Data is entered into the matrix row 
by row and fetched out columnwise. The inverse operation 
(deinterleaving) takes place in the receiver before the 
demodulation process. The details of the block interleaver 
implemented for this project are presented in Chapter II. 

Also important in mitigating the effects of fading 
for the satellite channel is finding the optimum modulation 
scheme in regard to bandwidth and power limitations. The 
typical satellite channel has the constraints of both 
bandwidth and power. One way to increase the power efficiency 
is to use error correcting coding by adding extra bits to the 


transmitted symbol sequence. This approach will decrease the 


bandwidth efficiency by forcing the modulator to operate at 
higher data rates which results in occupation of a larger 
bandwidth. For increasing the bandwidth efficiency, a higher 
order modulation scheme can be selected to maintain the same 
error probability; this will result in larger signal power 
consumption which is not power efficient [Ref. 4: p. 67]. 
Therefore, an optimum scheme must be found. 

Recent research on this problem emphasized the 
concept of combined coding-modulation schemes. In his 
important paper [Ref. 5], Ungerboeck states the advantages of 
Trellis coded modulation (TCM) as a tool that achieves 
significant coding gains over conventional uncoded multilevel 
modulation without compromising bandwidth efficiency. To 
implement this, one uses a convolutional encoder to increase 
power efficiency and finds a higher order signaling scheme to 
retain the bandwidth efficiency. For example, instead of using 
uncoded 4-PSK we could use a rate 2/3 convolutional encoder 
cascaded with an 8-PSK modulation scheme that will give the 
same throughput with some coding gain but without expanding 
the bandwidth required for 4-PSK. 

If we take into account the non-linear high power 
amplifiers which are generally used for satellite channels, we 
are led to consider chosing constant envelope modulation 
schemes. Continuous phase modulation has two attributes that 
make it attractive for the satellite channel: constant 


envelope and low spectral occupancy [Ref. 4:p. 240]. The 


application of coding to this modulation scheme adds memory to 
its own virtual memory to improve energy efficiency. This 
combined convolutional coding and continuous phase modulation 
Structure is of particular interest for energy constrained and 
band limited nonlinear satellite channels. 

Coded continuous phase modulation (CPM) was the 
starting point for this project. Existing satellite channels 
of interest use BDPSK modulation. We were to design a system 
that would have the same performance with a higher throughput 
or better performance with the same throughput. We chose a 
OPSK modulation scheme integrated with a conventional rate 1/2 
convolutional encoder. Both systems were tested through the 
fading satellite channel at different data rates. In addition 
we implemented minimum shift keying (MSK) modulation, a 
special form of CPFSK with modulation index h=1/2, through the 
satellite channel to compare its performance with those of the 
above two schemes. 

We Simulated the UHF Satellite channel by using the 
program PRO MATLAB on a SUN SPARK Station. The channel was 
under the effect of flat fading and AWGN. We used typical 
Space segment components like a prefilter, hardlimiter and 
postfilter to simulate the satellite transponder. 

Chapter II provides detailed specifications of the UHF 
satellite channel and its space segment components. Chapter 
III presents the simulation results for DBPSK, rate 1/2 


convolutional encoded QPSK and minimum shift keying (MSK) for 


different data rates. Detailed link budget calculations and 
conclusions are given in Chapters IV and V, respectively. 


Simulation codes are included as appendices 


II. UHF SATELLITE CHANNEL AND SPACE SEGMENTS 


A. FLAT FADING MODEL 


In this thesis we have used the Rayleigh fading channel 
model which is outlined in detail in Ref. 1, and shown in 


Fig. 1. 





Figure 1 Rayleigh Flat Fading Model 


The basis of the Rayleigh flat fading model, which is 
actually an implementation of a nonselective fading channel, 
is to impose a complex modulation in both amplitude and phase 
on the transmitted signal S,;=m(t)e™. Then we can represent our 


real received signal as 
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S,(t) =Re(A(t) e78'*) m(t) e7%*) (2-1) 


where A(t) and @(t) are time-varying amplitude and phase 
generated by the channel, m(t) is the transmitted data 
modulation and w, is the carrier angular frequency. This 
representation assumes that the fading spectrum is narrowband 
with respect to the carrier frequency [Ref.2: p. 52]. 
Furthermore, (2-1) can be represented in two quadrature 


modulating waveforms as follows 
Se) - omit) corw | +O. Cc) lm) cos (@ t+) ] BE) 


where 


DCG) =A ee) cost ( £) (4) 


Q(t) =A(t)sinO(¢c) 


These equations show that we can produce a flat fading 
received signal by passing the transmitted signal through a 
power splitter in case of non-complex working conditions, then 
shifting the phase of one of the outputs by 90 degrees and 
passing this output through two product modulators. The 
quadrature modulating waveforms I(t) and Q(t) are the channel 
disturbances. To simulate the channel accordingly, we have to 
define their statistics. 


As indicated in Ref. 1, for the Rayleigh case these I and 


alale 


Q modulating waveforms should be statistically independent and 
each should be zero-mean Gaussian process with equal 
variances. For our Simulation purposes we have used this 
model. For the DPSK case at baseband, the transmitted signal 
was real. We multiplied this transmitted signal with e” to 
have I and Q channels with equal power. For the other 
modulation schemes we worked in complex form so no power 
Splitter was needed. After we had our I and Q channels split 
from the transmitted signal, for each channel we produced 
white GausSian noise independent of each other and passed 
these Gaussian noise sequences through two cascaded RC filters 
whose specifications are defined in Ref. 2 as follows. 

Each RC filter had a time constant T,;=RC=7,/2.146193. These 
filters have been implemented by the simple difference 


equation 


Vy=AYVy.1 + DX, 


where y, is the output and x is the input of the filter. The 
filter coefficients a and b were determined by filter time 
constant and the filter gain respectively. When the filter 


power gain is set to unity, the values of a and b become 


= (2-4) 


PZ 


_ (1-a?)>, yi 
eer a (2-5) 


where A is the sampling time interval. These values of a and 
b affirm the resultant fading signal to have unit variance and 
zero mean. AS stated in Ref. 1, the value for A, should be less 
than or equal to 0.17, which is concluded experimentally. 
Therefore filter coefficient a has values between 0.2146193 
and 0.02146193. In our project we used 0.2146193 for fast 
fading and 0.02146193 for slow fading. 

The resulting fading signals are multiplied with the I 
and Q channel of the transmitted signal separately to produce 
the fading effected received signal. 


Figures 3 and 4 show the effect of slow and fast fading 


for a random generated input stream shown in Fig. 2. 


1 
o.5 | 
| 
—o.5 | 
| 

wn Q 300 1000 72300 2000 235300 pe Fo Loe] 335300 


time 





Figure 2 Transmitted Sequence 
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Figure 3 Slow Fading Effect 





Figure 4 Fast Fading Effect 


B. SPACE SEGMENT MODULES 

In addition to the fading simulation, we Simulated 
typical space segment modules for a satellite channel. The 
modules used during the simulation were adopted from Ref. 6. 
The modules were a prelimiter filter, a hard limiter and a 
postlimiter filter forv%a™ typical 25 Kizeenanne? } 

We will explain the main functions and detailed 


specifications of each below. 
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Le Prelimiter Filter 


The purpose of this filter is to establish the desired 
Channel bandwidth and provide an approximate 60 dB rejection 
for the stopband of 15 kHz from the corner frequency. The 
filter is implemented as Chebyshev Filter Type I with 6 poles 
and for 0.01 dB passband ripple. The magnitude and phase plots 
of the frequency response and step response of this filter are 


Shown in Figs. 5, 6 and 7, respectively. 
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Figure 5 Magnitude Plot of Prelimiter Filter 
Frequency Response 
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Figure 6 Phase Plot of Prelimiter Filter 
Frequency Response 
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Figure 7 Prelimiter Filter Step Response 


Filter function of the simulation program accepts the cutoff 
frequency value between 0 and 1, where 1 corresponds to half 
the sampling rate. For 12.5 kHz analog cutoff frequency the 


digital cutoff frequency was calculated to be 0.34 from 


2) a 


l curorF (normalized) a ‘r,N- z (2-6) 





Where Yo eeu. 
N : Number of discrete time samples per baud’ 
fae Analog sampling frequency 
2. Hardlimiter 
The hardlimiter is used to provide constant output 
power from the threshold of noise to maximum signal input. 


This was simulated by dividing each sample by its magnitude 


such that each complex sample is on the unit circle. 


$215) 


3. Postlimiter Filter 
The postlimiter filter attenuates out-of-band limiter 
noise power and spectral regrowth due to limiter and rejects 
multipath induced frequency response ripple. It was 
implemented as a Chebyshev Type I filter with 4 poles and 
0.025 db passband ripple. The magnitude and phase plots of the 
frequency response and the plot of the unit step response for 


this filter are shown in Figs. 8, 9 and 10, respectively. 
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Figure 8 Magnitude Plot of Frequency Response 
of The Postlimiter Filter 
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Figure 9 Phase Plot of Frequency Response of 
The Postlimiter Filter 
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Figure 10 Step Response of Postlimiter Filter 


Amplitude and phase responses of the channel are 
primarily determined by these filters. The low passband ripple 
requirement has imposed tight constraints over these filters; 
as a result, phase delays have occurred. An equalizer filter 
could be build to counterbalance this effect, but with 
Simulation tools we applied a certain amount of delay to the 
Sample stream to have the same approximate effect. Because of 
the hard limiter, the channel amplitude response is dominated 
by the postlimiter filter. Phase response, however, is 
dominated by both pre- and post-filters. 

With summaries of the flat fading model and the space 
segment units we concluded the presentation of the satellite 
channel. The overall UHF Satellite channel implementation 


diagram is shown in Fig. 11. 
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Figure 11 UHF Fading Satellite Channel 


C. MODULATION SCHEMES 

In this section the basic concepts of the modulation 
schemes that were tested with the UHF satellite fading channel 
are summarized. 

1. Differential Binary Phase Shift Keying (DBPSK) 

One of the important considerations in the modulation 
design is hardware simplicity. DBPSK modulation discipline 
does not require phase coherency and thus allows a relatively 
simple receiver structure due to absence of a synchronization 
Circuit [Ref. 7: p. 546]. This virtue makes it popular for 
channels where the phase shift changes slowly compared to the 
bit duration. 

Differential encoding is achieved by adding a 
reference bit at the beginning of the data stream. If the bit 


which comes after the reference is the same as the reference 


IL 


bit, no phase difference takes place. If the bit has a 
different value compared to the reference bit, 180 degrees 
phase difference applies; and so forth for each bit with 
respect to the preceding bit. 

There are primarily three decoding methods for the 
DBPSK as described in Ref. 8. We have simulated the bit-by-bit 
method of decoding DBPSK which theoretically requires, at 
most, 1 dB more Eb/No than that for Quadrature-PSK (QPSK) or 
BPSK provided that probability of error is 10% or less [Ref. 
73 Peo ocla 

Our implementation of DBPSK decoding was to take the 
difference and the sum of the current and previous symbol and 
compare the absolute values to determine the decision metric. 
Tf this metric is less than 0, then decision is made in favor 
of 0; otherwise 1. 

The modulator and demodulator implementation for this 


signaling scheme is shown in Fig. 12. 





Figure 12 DBPSK Modulator/Demodulator 
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We have taken this modulation technique as our base 
modulation scheme for the satellite channel simulation, since 
it is currently widely employed. Codes related with DBPSK are 
enclosed in Appendix A. We will present their performance in 
AWGN and fading channels in Chapter III. 

2. Rate 1/2 Convolutionally Encoded QPSK. 

We selected convolutionally encoded quadrature phase 
shift keying (QPSK) modulation scheme as an alternative for 
BDPSK, Since it has the same throughput and bandwidth as 
BDPSK. The convolutional encoder adds one more code bit for 
each bit in the data stream but, level 4 mapping of QPSK 
balances this excess bit so that the system has the same 
throughput and bandwith as DBPSK. Most important is that due 
to its convolutional encoded structure, this combined 
modulation scheme will be more immune than DBPSK to the 
destructive effects of the fading when blended with an 
interleaver. 

a. Rate 1/2 Convolutional Encoder 

Convolutional codes are named after the linear 
convolution that takes place between the information sequence 
and the code generator matrix to form the codewords [Ref. 9: 
me. 1). 

The encoder is linear and a finite state machine 


that has a sequential process ability. Therefore, it is 


AA: 


generally implemented by using shift registers. The number of 
memory cells in the shift registers plus one is equal to its 
constraint length K . The rate of the convolutional encoder is 
defined as the ratio of the number of the input bits divided 
by the corresponding output bits. 

In this project we have used the standard rate 1/2 
convolutional encoder which has (133,171), connections. The 
convolutional encoder and its connections are illustrated in 
Boao 

The maximum likely sequence estimator decoder for 
these codes uses the Viterbi Algorithm [Ref. 10]. The 
algorithm works iteratively symbol by symbol, tracing through 
the trellis to find the maximum likelihood correct path [Ref. 
8: p. 147]. The algorithm uses Euclidian distance branch 
metrics to compare each of the paths through the trellis for 


N symbol intervals (search depth) in the received vector. 





Figure 13 Rate 1/2 
Convolutional 
Encoder 
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b. Differential Quadrature Phase Shift Keying (DQPSK) 

OPSK modulation is a modulation scheme that uses 4 
phases per symbol. Thus, two bits are transmitted during each 
Signaling interval. 

QPSK signal is equivalent to two orthogonal BPSK 
Signals and it occupies the same bandwidth as that of a BPSK 
signal. 

For a satellite channel we desire a modulation 
scheme that does not require a synchronization circuit, nora 
noncohorent modulation scheme. For this reason we have encoded 
our QPSK signal differentially, as explained in DBPSK case. 
The performance of our implementation of DQPSK was 3 dB below 
the optimum matched filter coherent QPSK receiver [Ref. 7: p. 
548]. 

This combined convolutional encoded differential 
OPSK modulation scheme finally is blended with a block 
interleaver/deinterleaver to mitigate the effects of the 
fading channel by spreading out the error bursts in time 

The codes related to OPSK modulation simulation are 
in Appendix B. 

c. Block Interleaver 

The block interleaver/deinterleaver for this 
application was basically a two-dimensional matrix with k rows 
and 1 columuns where k.l=Total number of symbols. The number 


of rows (k) determines the degree of the interleaver and is 


Zs 


set long enough for slow fading to ensure the errors within a 
code word are independent. However, the delay caused by this 
remapping in real modem implementations must be taken into 
consideration. 

The block diagram of the overall implementation for 


the combined structure 1S Shown in Fig. 14. 





Figure 14 Block Diagram of the Rate 1/2 Encoded 
DOPSK Modulation Structure 


3. Minimum Shift Keying (MSK) 
As we stated in Chapter I, there is an interest in 
constant envelope modulation schemes in the satellite channels 


due to nonlinear high power amplifiers. Continuous-phase 
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shift keying (CPFSK) modulation has a constant envelope 
Structure and restricted spectral width. We simulated this 
modulation scheme and passed it through fading channel to have 
a brief understanding of its spectral efficiency in the 
satellite channel and its response under fading conditions. 

The continuous-phase shift keying signal is generated 
by feeding the data signal into a frequency modulator. The 
output of the frequency modulator will have signal shape 


described by [Ref. 3: p. 173]. 


s(t) =Acos[2nf.t+®(t;I)+,] (227) 
where 
A : Signal amplitude, 
in : Carrier frequency, 


@(t;I) : Time-varying phase of the signal, 


do : Initial phase, 
fe : Time, 
af : Modulation value. 


The time-varying phase of the signal can be described as 


oak 
O(f; I) -2nt,7 d I,+2nf,q(t-nT)I,;nT<ts(n+1)T (2-8) 


or by equating h(modulation index)=2f,T we can rewrite the 


above equation as 


Agt 
b(t;I)=nh XZ I,+2nf,1,q(t-ntT) (259) 
kK=-© 


ZS 


where 


ee : Peak frequency deviation, 
T “Symbol durarven, 

n : Symbol number, 

ia : Modulation value, 

eileen : Phase shaping waveform. 


The phase shaping waveform is the integral of the frequency 


Shiite, (Chatweus, 


C 


a(t)=fu(t) de (2-10) 


0 


The frequency pulse u(t) has a finite duration for the 
interval OstsLT. L is the pulse length and determines the 
response of the modulation. L=1 is called full response CPFSK, 
and L>1 is partial response CPFSK. Although frequency pulses 
have discontinuities, the integration makes the phase-shaping 
pulse continuous, and the resulting waveform has a continuous 
phase which is a function of both time and modulation values. 
In (2-8) I=....1,,1,,1),1;,..-. denotes the modulation sequence 
sent to the CPFSK modulator. I, takes values {+1,+3,....,+(M- 
1)} where M =2 for M-ary CPFSK, and k corresponds to number 
of information bits per symbol. The parameter h (modulation 


index) determines the number of states in CPFSK. If his a 
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rational number, then the number of states in CPFSK is finite. 


It 1S common practice to select a modulation index 


(2 =a) 


= 
iQ Io 


where p and g are relatively prime integers. These integers 
determine the phase state so that when p is even, there are g 
Sates and 2q states when p is odd [Ref.3: pp. 172-175]. 
CPFSK signals are generally characterized by their phase 
trajectory versus time. These phase trajectories can be 
represented in a trellis structure, which for the binary 


sequence is shown in Fig. 15. 





Figure 15 Phase Trellis of Binary CPFSK 
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Changing three parameters of the CPFSK' structure, 
modulation index, pulse shape and the alphabet size of the bit 
stream, enables one to generate variety of CPFSK signals. One 
important of these is MSK signaling, which is binary CPFSK 
modulation with index h=1/2. MSK modulation uses minimum 
frequency separation of Af=1/2T to ensure that the signaling 
waveforms are orthogonal [Ref. 3: pp. 180]. 

In our simulation, the MSK waveform was produced by 
implementing (2-6) for baseband signal generation with zero 
carrier frequency, as described in [Ref. 7: pp. 351-353]. The 
demodulator used for the simulation was an implementation of 
the FSK demodulator for an MSK waveform [Ref. 8: pp. 258-259]. 
In the demodulator we used two match filters to extract the 
binary signals and applied their output to the Viterbi 
algorithm for soft decoding. Appendix C contains simulation 


programs related to MSK modulation. 
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III. SIMULATION RESULTS 
In this chapter we present the simulation results for 
DBPSK, DQPSK and MSK modulations in the fading satellite 


channel. 


A. DBPSK SIMULATION RESULTS 
The simulation for DBPSK modulation was conducted in three 
categories. 
1. AWGN Channel Simulation 
We used 10000 bits to simulate the base DBPSK 
modulation and compared this with the theoretical curve to 
make sure our Simulation setup was accurate. The results are 


illustrated in Fig. 16. 





Figure 16 Theory and Simulation 
Comparison for AWGN Channel 
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The simulation results agree with the theoretical 
results. The theoretical curve was plotted using [Ref. 7/2gee 
549] 

522) 
Pi=se No (3 =a5) 


e 


We plotted our simulation points as the ratio of number of 
decoding errors out of the total number of transmitted data 


bits versus E,/N, , where 


n=N 
2 
(enh ain Set (3-2) 
‘N, 2k g2 


and A is the signal amplitude, N is the number of samples 
Guring a symbol period, k is the number of bits per symbol and 
o is the standard deviation of the AWGN. In our implementation 
the values chosen were Ae {1,-1}, N=30,k=1. The o varied 
according to the E,/N, under test. 
2. AWGN Channel and Transponder Components 

In this test, we have included the prelimiter filter, 
hard limiter and postlimiter filter portions of the satellite 
transponder to determine their effect on the performance. We 


conducted the simulation for 2400 baud with the satellite 
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transponder components included. The results are shown in Fig. 


a . 





Figure 17 Effect of Transponder Components 


The data show that the effect of inclusion of 
transponder components in the channel is approximately 0.8 dB 
loss in performance. We used the same setup to see the 
response of 4800, 9600 and 19200 baud data rates by changing 
the number of samples per symbol and corresponding noise 
Standard deviation. These data are shown in Fig. 18. There is 


almost constant 1 dB loss for each data rate increase. 


fal 





Figure 18 Different Baud Rates for DBPSK cascaded 
with Transponder Components 


3. Satelite Fading Channel 


We used the fading generator in the slow fading mode 
as explained in Chapter II for this simulation. We tested four 
transmission rates through this channel as shown in Fig. 19. 

If we compare the theoretical performance for AWGN 
channel with the 2400 baud fading effected performance we see 
a catastrophic degradation which requires higher signal power 
for a desired bit error rate. This is an unacceptable solution 
for satellite channels, due to physical constraints of the 


space segment. 
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Eb/No (dB) 





Figure 19 Fading Satellite Channel Effect over Different 
Baud rates for DBPSK Modulation. 

B. DQPSK SIMULATIONS 

Considering the simulation results for the fading 
Satellite channel and noting that the energy requirement is 
very high for desired bit error rate, we proceeded to 
implement rate 1/2 convolutional encoded DQPSK modulation for 
the channel. Simulations related to DQPSK were performed in 


the following categories. 


33 


1. DQPSK Simulation in AWGN Channel 
To be sure that our DQPSK simulation codes were 
working properly, we checked DQPSK in AWGN channel. The 
outcome of this simulation is plotted in Fig. 20. For this 
Simulation we have used k=2 in (3-2) for the relation between 
E,/N, and the standard deviation (0) of the AWGN. The 
theoretical curve was plotted 3 dB above the performance of 


coherent QPSK, which is 


4p (Gey 
P.=0(,)2(=2)) 


.o) 


The figure shows that simulation results matched well with the 


theoretical curve. 





Figure 20 DOPSK Simulation in AWGN 
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2. Rate 1/2 Convolutionally Encoded DQPSK in AWGN Channel 
This second category simulation was also a test 
Simulation to check the convolutional encoder performance. In 
this case we used k=1 in (3-2) for the bit energy per 
information bit transferred. The signal constellation at the 
output of the differential demodulator for 13 dB E,/N, is 


illustrated in Fig. 21. 





Figure 21 Signal Constellation of Convolutional 
Encoded DQPSK in AWGN Channel for 13 
dB E,/N,. 


ie, 


The simulation® reswits “for the convolutioename, 
encoded DQPSK in the AWGN channel are shown in Fig. 22. In the 
figure, BER slowly changes until it reaches the threshold of 
10°, then the effect of convolutional encoding becomes 
dominant; we see an asymptotic performance curve as is 
characteristic of systems with forward error correction 


Codimg: 





Figure 22 Rate 1/2 Encoded DQPSK in AWGN Satellite 
Channel 
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3. Rate 1/2 Convolutionally Encoded DQPSK in Fading 
Satellite Channel 

Given the results in AWGN channel we, next cascaded 

Our transponder components and fading generating function to 

the simulation program and tested the entire system for 

different baud rates. Figure 23 presents the signal 

constellation at the output of the differential demodulator 

Bereetiis system. Comparison of Fig. 23 and Fig. 21, 

constellation in AWGN, shows the effects of fading are 

principally on the signal amplitude which now varies in 


accordance with the Rayleigh probability distribution. 





Figure 23 Signal Constellation for DQPSK 
in Fading Satellite Channel for 25 dB E,/N,. 
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The constellation iS »£OnsEp/Ne2S.dB. Thoughethe err-—s 
of the noise is very small, the amplitude variation due to 
fading is coarse. This variation is the source of the 
performance degradation due to fading. 

The performance of the convolutionally encoded DQPSK 


is illustrated in Fig. 24. 





Figure 24 Performance of Rate 1/2 Convolutionally Encoded 
DOPSK in Fading Satellite Channel 
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This simulation was done with transponder components, 
including interleaver/deinterleaver. Interleaving depth is 
0.01 of the total transferred bits. Fading was applied to make 
differential phase shift between symbols on the order of 10° 
radians. We plotted the curves related to baud rates 2400 and 
19200, due to very close simulation outcomes. 

If we compare these performance curves with those of 
the AWGN channel we see an approximate performance 
degradation of 13 dB. This loss is due to fading effect of the 
channel and is the performance that can be achieved with 
present fading mitigation methods. But this performance, for 
the same throughput, is superior to the performance of DBPSK 
in the fading channel which is entirely unacceptable. (See 


Fig. 19.) 


C. MSK SIMULATIONS 

We simulated MSK modulation for the satellite channel in 
two categories: first, for the AWGN channel with transponder 
components included, and then with the fading portion of the 
code included to observe the effect of fading on MSK. 

1. MSK Performance in AWGN Channel 

We used 10000 bits in these simulations. The E,/N, was 

calculated as described in (3-2). The theoretical curve of 


the MSK in AWGN channel is [Ref. 7: p. 549], 
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i Eb (3-4) 
as ol 7am 


The performance shown in Fig. 25 is the result of simulation 


in the AWGN channel with transponder components included. 





Figure 25 MSK Performance in AWGN Satellite Channel 
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2. MSK Simulations in Fading Satellite Channel 
This second type of simulation was done in the 
presence of fading with transponder components included. The 
performance of MSK in the fading channel was severely degraded 


compared to that in the AWGN channel. 


Theory in AWGN 
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Figure 26 MSK Performance in Fading Satellite 
Channel 
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IV. LINK BUDGET MARGINS DERIVED FROM THE SIMULATIONS 
In this chapter we will show the link budget margins 
between the calculated values of the UHF Follow-On Satellite 


System and the simulation results. 


A. UHF FOLLOW-ON SATELLITE SYSTEM CALCULATED VALUES. 
The parameters of calculating the carrier-to-noise ratios 


are defined in Ref. 6. These values are tabulated in Table 1. 


TABLE 1 


UHF FOLLOW-ON SATELLITE SPACE SEGMENT SPECIFICS 





27 dBW 29 dBW “14.4 dB/K 


where 
EIRP = Effective Isotropic Radiated Power, 
G/T = Receiver Sensitivity. 


For evaluation purposes we chose the AN/WSC-3, the US 
Navy’s primary mobile UHF satellite terminal transceiver whose 
characteristics are defined in [Ref. 11: p. 16] and shown in 


Table 2. 
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TABLE 2 


AN/WSC-3 TRANSCEIVER SPECIFICS 


22 dBW Sle St toe 1s 


We selected two extreme cases for the location of the 






terrestrial terminals. The first terminal was under direct 
projection of the satellite coverage which has an elevation 
angle of 90 degrees, and the other was with the minimum 
elevation angle of 5 degrees. For these stations we 


calculated their slant distances as shown in Table 3. 


TABLE 3 


TERMINAL SLANT RANGES 


eee ELEVATION (E) SLANT DISTANCE (d) 


357/68 sn: 
ee. = 


The corresponding free space path losses are calculated as 











Er ae 
| te 





follows; 


2Fy0h og as (4-1) 


é = 
(dB) Cc 


EL 
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Where 


f 


d 


Cc 


¢c 


Carrier frequency (For our case 300 MHz.) 


Slant distance (km) 


Speed of light (2.99 x 10° km) 


The free space path losses for stations 1 and 2 were 


calculated and found to be 173 dB and 174.29 dB, respectively. 


The carrier-to-noise ratios for uplink and downlink cases 


were calculated by using below equations [Ref. 11: p. 74]. 


Where 


a 


one, 


N 


oS 
F 


Bian Dieea 
EIRP seme 


G/ a eniinsl 
G/ Bt rareliite 


BW 


G 
) AOR eae come ( a terminal 


G 
=! uptaiieee 1 tee eee ( = satelli cael O ec cai = Bie (4 -2) 


-10log,,k-BW (4538 


EIRP of the ground terminal (dBW), 

EIRP of the satellite (dBW), 

Receiver sensitivity of the ground 
terminal (dB/K), 

Receiver sensitivity of the satellite 
(dB/K), 

Noise Bandwith (1010g,)(25 kHz)), 

Free space path loss (dB), 


Boltzmann's constant (dBW/K-Hz). 


ab 


Calculated carrier-to-noise ratios’ for uplink and 
downlink are tabulated in Tables 4 and 5 respectively. 
TABLE 4 


UPLINK CARRIER-TO-NOISE RATIOS 


STATION (ON) opin 






TABLE 5 


DOWNLINK CARRIER-TO-NOISE RATIOS 





B. SIMULATED CARRIER-TO-NOISE RATIOS 

The performance curves of the simulations performed in the 
Satellite channel gave us the E,/N, values for corresponding 
probablity of error. By using (4-4) we converted our §£,/N, 


values to carrier-to-noise ratios for different data rates. 


a) 


Sle 
| 
by | by 
Id 


45 


Where 
R = Baud rate (symbols per second), 


B Noise Bandwith (25 kHz). 


From this relation we calculated our simulated required 
carrier-to-noise ratios for BER 10* and tabulated them in 
Table T6e 

TABLE 6 


SIMULATED CARRIER-TO-NOISE RATIOS 


AWGN -0.17 dBwl 363° aBe 4 7 oee~ ane 118s 











For MSK in the fading satellite channel we did not 


calculate the corresponding carrier-to-noise ratios, because 
they were out of range for BER 10%. 

We calculated our carrier-to-noise ratios for 19200 baud 
rate, but it should be noted that maximum baud rate for 


AN/WSC-3 is 9600 baud. 
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C. LINK MARGINS 

If we compare the contents of Table 6 and with those of 
Tables 4 and 5 for the fading satellite channel, we see that 
the only suitable modulation technique is convolutionally 
encoded DOPSK, which is lower than the carrier-to-noise ratios 
for station 1 in all cases and for station 2 in all cases 
except baud rate 19200 in the relay mode of the satellite. 

We tabulated link margins in Table 7 for the enhanced mode 
of the satellite and for the worst case of earth station, 


Station 2, which has more slant range. 


TABLE 7 


LINK MARGINS FOR ENHANCED MODE 


MODULATION 2400 BAUD 4800 BAUD 9600 BAUD | 19200 BAUD 


| DBPSK AWGN Dimoge ar ti gmome de l2e7600dR 1) o=25 fap 
DBPSK FADING oy GBM eet Om 7 ede) - 17 746d - 22°75 "dB 


AWGN gy Al O 4 wert 7 ds Peo coo Deas cs 





For the other fading situations, predicted carrier-to 
noise ratios are greater than the actual proposed carrier-to- 


noise ratios. Therefore, there 1s no margin for these cases. 


4'7 


For the AWGN channel, in all cases the simulated carrier- 


to-noise ratios are below the actual values of the proposal. 
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V. RESULTS AND RECOMMENDATIONS 


For this thesis we simulated the satellite fading channel 
proposed in Ref. 6. We also simulated DBPSK, rate 1/2 
convolutionally encoded DQPSK and MSK for both AWGN and fading 
satellite channels. We evaluated their performances for these 
channels for different baud rates and for different BER 
values. 

The results show that the fading effect of the channel 
degrades the performance of the DBPSK severely. We tested a 
convolutional encoder and interleaver/deinterleaver pair 
cascaded with DQPSK modulation to mitigate the effect of 
fading without changing the throughput relative to DBPSK. We 
see that the performance of the encoder combined with DOQPSK is 
20 dB better than that of the DBPSK in the fading satellite 
channel. Additionally we simulated MSK for different baud 
rates in AWGN and saw that it has an approximate 1 dB loss due 
to baud rate increase. However, its performance in fading was 
severely degraded. 

Our ultimate aim was to cascade this MSK modulation with 
a convolutional encoder and apply some diversity technique to 
evaluate its performance for the satellite channel due to its 
smooth spectral character. Follow-on studies may provide 


opportunities for other students to combine encoding and 
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continuous phase modulation (CPM) with the satellite fading 
channel. However, how to interleave coded CPFSK symbols is an 
unsolved problem. 

Fading also can be modified to simulate the multi-path 
effect for present modulation techniques and combined coding 
and CPM. 

Further, two channels should be simulated simultaneously 
to determine the effect of adjacent channel interference on 


the different modulation and coding schemes. 
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APPENDIX A. DBPSK CODES 


A. DBPSK MAIN PROGRAM 


elear; 
omg 


Sok KKK KKK KKK KKK KKK KKK KEKE KK KEK KKKEKEKKEKEKKKEKKEK 


Definition of variables 
SR KKK KEKE KEKE KE KEK KKK KKK KKK KKK KKK KEK KEKE KEKKEKEKKEKKEK 


m=10000; sNumber of bits 

del=.01 sDifferential phase shift for the fading 
or (30 15 7 4); «Number of samples per symbol 

Sagia= 1.5 ; *Standart deviation of the noise 


filt coeff pre=.32; *Prelimiter filter coefficient 
fFilt coeff post=.21;7Postlimiter filtSer coefficient 


RakkkKkKkKKKKkKkKkKKKKKKKKKKKKKKKKKKKKKKKKKKEKKKE 


For each sample per symbol and for each noise standart 


deviation program runs once and saves the number of errors 
SKK KKK KKK KKK KEKE KKK KEKE KEKEKKEKKKKEKKKKKKKKEKEK 


for kk=1:length (cf) ; 
for i=1:length(sigma) ; 


feet (KK); 

n coef=sigma (1) ; 

rnd_o=msg(m) ; *Random bit stream creation 

See o—Gitecod(rnd 0); tDifferential decoding 

map o= map(dif_o); Map {0,1} to {-1,1} respectively 


mod _o=modul(map_o,N); *Phase Shift Keying Modulating 
mod _o=mod_o’; 
meq O=mod 0(:) ; 


mod _o=exp(j*pi/4) *mod_o; sPower Splitting 

I=fading (del, (m+1) *N) ; $I channel fading stream creation 
Q=fading (del, (m+1) *N) ; *Q channel fading stream creation 
mod_o=I.*real(mod_o) +j*(Q.*imag(mod_o)); *Fading addition 


save ml mod_o; 

clear rnd_o dif_o map_o mod_o I Q; 

ch_dis=awgn_ch((m+1)*N,n_coef);sAWGN creation 

load ml; 

dem_in=ch_dis+mod_o;sAWGN addtion to the modulator output 

clear mod_o ch_dis ; 

buf=ones(2*N,1) ; 

dem_in=[buf;dem_in];%Filter steady state time compensation 

(b1,a1] =cheby1(6,.01,f1i1lt coeff pre) ;%Prefilter 
scoefficients 


Sb. 


dem in=filter(bijjal,dem in); “sPreti leew 


dem_in=limiter(dem_in); % Hardlimiter application 
[b2,a2] =che@by1 (4, .025, fit ceoeit pest) -Pest caaeer 


coefficients 
dem in=filter(b2,a2,dem_in);%*Postfiltering 
dem_in=dem_in (2*N+1;length(demiin)), :);sshitting ducmae 
filter 


dem_in=[(dem_in(8:length(dem_in)) ;dem_in(1:7)] ;%delay 
dem in=reshape (dem iin, Nine); 

dem_in=dem_in’; 

[out, metric] =demod(dem_in,N,m) ;’Demodulating 

clear dem_in; 

save ol out; 


load ri; 
[n_er(kk,1)] =compa(rnd_o,out,m) ;sComparison of received 
sand transmitted sequences 
end 
end 


Save res.1 n_er %Saving number of errors 


B. MSG.M (INPUT SEQUENCE GENERATION FUNCTION) 


function u = msg(k) 

This M-file accepts k the number of bits that will be 
returned in the vector u for seed 1. 

rand (’ unttonmns 

Wo = LOundi rand (ioe, 


C. DIF COD.M (DIFFERENTIAL ENCODER FUNCTION) 


functions (dit ol=dit (coda) 
This M file accepts an input vector and differentially 
sencodes it by inserting the reference bit 1 at the 
sbeginning of the vector. 
a=length(in) ; 
y=[(1,zeros(1,a)]; first 1 is the reference bit 
for i=l:a; 
(tee) ) = en ta eye eee 
end 
dif o=y; 


D. MAP.M (MAPPING FUNCTION) 


Function o=map (in) 

eThis M file accepts the binary vector and mappes 0’s to 
‘1's. 

a=length (in) ; 

fous? —ilea 
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end 
O=in; 


E. MODUL.M (MODULATOR FUNCTION) 


mde: lon [y] =modul (in,N) 

«This M file accepts the differentially encoded and mapped 
evector in and N number of samples per symbol and and 
sgenerates a matrix with dimensions (length(in) x N). 
a=length (in) ; 

y=in’ *ones(1,N) ; 


F. FADING.M (FADING GENERATION FUNCTION) 


function y=fading (del, mm) 

6This M file accepts differential phase shift del and total 
sénumber of the samples mm in the input vector and creates 
4mm normal distributed variables and passed these through 
6two RC %filters. 

s=exp (-del/2.146193) ; 

ss=((1-s*2)*3/(1+s°2))*.25; 

b=[ss]; 

a=[1 -s]; 

rand(’normal’ ); 

G=rand (mm, 1); 

k=filter(b,a,t); 

metilter(b,a,k); 


G. DEMOD.M (DEMODULATOR FUNCTION) 


function [out,a] =demod (mod_o,N,m) 
¢This M file accepts fadign and AWGN effected modulator 
‘Output mod out, number of sample per symbol,and total 
snumber of symbols. And for each symbol takes the difference 
sand sum of the present and and past symbol sums and takes 
squares of their absolute values.And if dif < sum then 
6decides input bit is 0. 
o=ones (1:m); 
mor L=2>m+1: 
dif (i-1)=abs (sum(mod_o(i,:))-sum(mod_o(i-1,:))).%2; 
su(i-1)=abs (sum(mod_o(i,:))+sum(mod_o(i-1,:))).%2; 
metric (i-1)=dif (1-1)-su(i-1); 
1f metric(i-1) <0; 


3 


O(1-1) =0; 
end 
end 
a=metric; 
out=0; 





APPENDIX B. DQPSK CODES 


A. DQPSK MAIN PROGRAM 


clear 

elg 

@el—0.1; 

m= !O000; s¢Number of symbols that wil be transmitted 
eaeaEKKKKKKKKEKKKKKKKKKKKKKKKK KKK KKKKKKKKKKEKKKKKKEKE 

6You can change number of samples N accoeding to the 
sBaud rate you desire to work on.N=30 for 2400,N=15 
for 4800,N=7 for 9600 and N=4 for 19200 Baud.You have 
6to consider to change the standart deviation of AWGN for 
seach case. 
KKEKKKKEKEKEKKKKKEKEKEKEKKKKKKKKK KK KKKEKKKKK KKK KKK KE KKEKEKEKSE 


N=30;%Samples per symbol (This is for 2400 baud) 


load fec Data file that has G64 and Té6é4 

u=msg(0,m); *Random symbol creation 

sesoxrt (-1); 

Save utempl u 
v=cnced(2,1,6,G64,u, zeros (1,12)) ;sConvolutional encoding 


mm=bm(2,v); ¢Binary to M-ary conversion 
clear v u; 
mm=inter(50,200,mm) ;%Interleaving (50 by 200) 
M=dphmod(2,mm) ;#Differential encoding 
clear mm ; 
M=resamp(N,M) ; Generate each element of M , N times 
M=reshape (M,N,m+1);%Reshape to build (m=1) by N matrix 
M=con] (M’ ); 
save Mtemp M; 
clear M; 
sigma v=[{1.2247 ];%sAWGN standart deviation 
for ss=1:length(sigma_v) 
Np=80; % This is the path history length for the 
6viterbi decoder 
KK=fad(del,m+1,N) ;%Applying fading 
load Mtemp; 
M=KK. *M; 
clear KK; 
R=awgn(M,sigma_v(ss));sAWGN addition 
Re @@lapia }; 
R=R(:);¢Make it vector to pass it from filter 
clear M; 
buf=[{ones (2*N,1)]; 
R= VouUreR| - 
(ede osie = Cheby AGO | 32). 
RS amewer (old l, ke); 


oe) 


R=limiter(R); 
[b2,a2}=cheby (4, 02572 
R=filter(bZ,agan: 
R=R(2*N+1:length(R)) ; 
Ree 7 


R=(R(@;lengcen (Rien, 

R=reshape(R,N,m+1); 

R=conj (R’ ) ; 

for 1=1:m+l 
temp (i, :) =xcorr(ones(1,N) *exp(j*pi/4),R(i,:)); 
6This correlates each symbol with ones to find 
metrics 
res (i) =(temp(i,N))/N;%Normalization 

end 


clear temp R ; 
[s,R] =ddmod(2,2,res) ;%Differential demodulator 
clear res Ss; 
R=dinter(50,200,R); *Deinterleaving 
D=eucdis(2,R);%Euclidian distance finder 
clear R temp2; 
[rr,cc] =size(D) ; 
PH=zeros (64,3*Np) ; 
fOr, a lates 
PH=viterbi(1,Np,PH,T64,D(h,:));%Viterbi decoding 
ud (h)=PH(1,240)-1;%Signal level offset 
end 
clear PH D;load utemp; 
[a,b] =mtrxv(79,u,ud);%Shifting due to path history 
$length. 

clear u ud; 
err(ss)=check(a,b);*Error check 
clear ab; 

end 

Save result err 


B. CNCD.M (CONVOLUTIONAL ENCODER FUNCTION) 
function ([v,VtCle= coc no. ne Ge ee 
CONVOLUTIONAL ENCODER 
Paul H. Moose 
Naval Postgraduate School 
04-24-91 


CP olf off off ol? of? 


This m-file is a feedforward convolutional encoder for 
san(n,k,m) linear convolutional code.It codes binary vector 
4u of length Lk bits into binary vector v of length Ln bits 
¢plus a reSidual vector vr of nm bits.The code is determined 
sby the matrix filter Gr of dimension k rows by n(m+1) 
¢columns. (See page 292 in Lin and Costello Gr=[{G0O Gl G2 
sGm].). 

S A previous initial vector r of nm bits is added to the 
6first nm bits of v. 


9] 


a6 


The program first checks to see if Gr has correct 
Semeislons el: @t goes, 1f continues. If net, itwends and 
sreturns an empty v. 

ikg N)]=size(Gr) ; 


of o\? 


1£ N == n*(m+1) 
Matrix has correct no of columns 
oy == 


Matrix has correct no. of rows 
«The program now checks to see if u has a multiple of k 
soits.if so, it continues, if not it adds zeros to u to make 
Sit so. 
[zz Lk] =size(u) ; 
gq=rem(Lk,k) ; 

im == 

u=uU; 

else 

u= [u zeros(1,k-q)]; 

[zz Lk] =size(u); 

end 
% We now have a msgword u with Lk bits, an exact multiple of 
$k bits. The output code word is formed by finding first the 
‘output for each word of k bits. These outputs are make the 
Srows of w. 
L=round (Lk/k) ; 


Bem j=1:4 
wae, 2) =rem(u(:, (j-1) *k+1:]*k) *Gr, 2); 
end 


The output vector is formed by adding rows of w with 
sprogressive delay of n bits for each row. The initial 
condition vector r is placed in the first nm bits of v. 
v=(r zeros(1,n*L)]; 


ote =e 1) 
v(:,n*(1-1)4+1: (m+1) *n)= v(:,n* (1-1) +41: (m+1)*n) +w(1,:); 
end 


v=rem(v,2) ; 

¢This vector is divided into the residal portion of nm bits 
sand the first nL bits. 

vr=v(:,n*L+1:n* (L+m) ); 

amt; , 1: n*L) ; 

6This is the end of the program for v. 

else,break,end %End if wrong no. of rows 

else,break,end %End if wrong no. of columns 


C. BM.M (BINARY TO M-ARY FUNCTION) 


PaAGeLOn m = bm(q, v) 
BINARY TO M-ARY CONVERTER 
Paul H. Moose 
Naval Postgraduate School 
SS) = ALO Shab 


a oP olf ol? Fh 


oa) 


The bits are stripped q at a time and are mapped 
to a symbol vector m with integer values 0 to 2%q-1. 
zeros are added at the end of vy ,to have aGmgeqeq = 


a 

o 

& 

[>] 

a 

fo 

46 if necessary. 
N=Z ea. 

[m n] =size(v) ; 
1G 
7 


= rem(n,q) ; 
f Yr == 
V=V; 
else 
Vv = [v Zercst® g-=) |; 
end 


* We now have a codeword v with 
¢ The bits must now be stripped 
s6generate the m-ary values. 

map = 1; 

for j =l:q-1 
map =[map 2°4]; % 

* the left 
end 

[m n] = size(v); 

Dp) =neunedin, oF; 

her r= ita 


s = v(:, (1-1) *q+1l: (1-1) *q+q); 
m(i) =s*map’; 
end 


an exact multiple of q bits. 
gq at a time and used to 


This makes the least significant bit on 


D. INTER.M (INTERLEAVER FUNCTION) 


PtactlonsmaiH1n tava mh) 


*This m-file is a block interleaver. The matrix m is read 
4into a column by rows, then the column is read into 


6 a (1,kK) matrix by columns and 


read out into vector mi by 


SrowS. mi is shaped into matrix of same dimension as Mm. 


M=—Gomyem.) 7 
y=m(:) ; 
N=length(y) ; 
eee ny 
= Zens: le hae 
x(:) =y; 
X=COngoe 
VE=7e ee 
m(:)=yi; 

Ma —COma ins) 
end 


5 )5) 


E.DPHMOD.M (DIFFERENTIAL PHASE MODULATOR) 


pom@etaon MD = dphmod(q,m) 

sThis M-file creates complex values MD with amplitude one 
stand one of 2°q equal phase values. The first value is one. 
eThe next values are differentially coded in phase. The 
Sinput symbols are in the m vector. 

N=2°q; 

Gen=2*pi/N; 

[rr n]=size(m); 

Ber k=l1:rr 


md=cumsum(m(k,:)); This differentially codes the phase 
values 

e—Sort (-1); 

MD(k,:) = exp(i*dph.*md) ; 

end 


MD=[ones(rr,1) MD]; 


F. RESAMP.M (SAMPLING FUNCTION) 


function y =resamp (k,x) 
This M-file resamples the vector x to generate y. If k is 
snegative,every kth sample is put into the vector y. If k is 
positive, every sample in x is repeated k times in y. 
[m n] =size(x); 
if k >0 
mom, 1=):n 
Bem i= 1 :k 
waa * i - 7 +11) =x (age: 
end 
end 
else 


G. FAD.M (FADING FUNCTION) 


function y=fading(del,aa,bb) ; 

This M file accepts del, the differential phase shift and 
6the dimensions of the modulator output and calculates 
sfilter coefficients for 2 RC filters, and generates 2 
snormal random distributed sequences passes these sequences 
éthrough the filters and forms the Rayleigh fading signal. 
s=exp (-del/2.146193) ; 
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ss=((1-s°2)°3/(1+s 2 
b=(ss]; 

a=(1 -s]; 
rand(’normal’ ) 
rand('seed’ ,0O); 
k=rand(aa*bb,1) ; 
k=filter(b,a,k,1); 
k=filter(b,a,k,1); 
z=rand (aa*bb,1) ; 
Z=feel ter( bp a;Z el): 
z=filter(b,a,z,1): 
k=reshape (k,aa,bb) ; 
Z= ee? eae aa, bb) 
jJ=sqrt(-1); 
y=k+j*Z; 


i 


H. AWGN.m (AWGN FUNCTION) 


function y = awgn(X,sigma) 

* Awgn is an m-file that adds awgn noise to the matrix x 
% Eb/No=1/(2*sigma~2) . 

[rr , cc] =e Zevear, 

rand (’normal’ ) 

w=rand (rr,cc)+i*trand(rr,cc) ; 

y=x+Sigma. *w; 


I. DDMOD.M (DIFFERENTIAL DEMODULATOR FUNCTION) 


function [s@M) =ddmoed(qp, a, MD) 

This M-file differentially demodulates complex modulation 
svalues in MD into 2°qp equal phase sectors from 
constellations of 2°q phase sectors 

ssectors into ka amplitude bits. the output is [s M] where 
6 S is the phase sector number and M is the differentially 
+decoded modulation values. 

k=length (MD) ; 

£Om j—iec- 1 

(a)) =MB (+2) * cong wigs, 

aan 

N=2°qp; 

aon ah 

phn=angle(M) ./dph; 

s=rem(round (phn) +N,N) ; 


S 


60 


J. DINTLV.M (DEINTERLEAVER FUNCTION) 


function m=dintlv(1,k,mi) 

¢6This m-file is a block deinterleaver. The vector mi is 
sread into a (k,l) matrix by columns and read out into 
evector m by rows. 

me —cony (mi’) ; 

y=mi(:); 

N=length(y) ; 

ees k==N 

x=zeros(k,1); 

x(:)=y; 

mecOny (x'); 

bel=x(:); 

aie: ) yi; 

m—con) (mi’ ) ; 

end 


K. EUCDIS.M (EUCLIDIAN DISTANCE FUNCTION) 


humet ion Di = eucdis(q,R) 
EUCLIDEAN DISTANCE METRICS 
Paul H. Moose 
Naval Postgraduate School 
06-17-92 


This M-file finds Euclidean distance of elements 

in vector R from 2°q unit amplitude vectors 

equally spaced on the unit circle. It stores these as rows 
eOE D. 

N=2°q; 

f= hength(R) ; 

index=1:N; 


HP ol? ol off olf ol? Ol? ol? 


Oo 


siem—<*pi/N; 

MO=exp (j *dph.* (index-1)) ; 

mom l=): 
D(1,:)=abs(R(1).*ones (MO) -MO); 
end 


L. VITERBI.M (SOFT VITERBI DECODING FUNCTION) 


function PHN = viterbi(k,Np, PH,T,D) 
SOFT VITERBI DECODER 
Paul H. Moose 
Naval Postgraduate School 
LOC 


This m-file decodes k bit msgwords from 2*n real 
metrics in the row vector D. (These may, for example, 


a? o\P olf olf ol9 olf ol? 
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srepresent the "Euclidean distance" of the received 
modulation value from each of 2°n modulation values.) 

The state transition information for a K state trellis is 
stored in a K by 3*2°k matrix T. Each of the 2*k paths 
sentering each state has its source state (one of K), path 
smsgword (one of 2°k) and path codeword or signal number 
(one of 2*%n) listed in the state row of T. 

S The Np length path histories are kept in matrix PH which 
is K by 3*Np. The path history for each state contains 
ssource state, path weight and path msg word for Np 
sprevious states. 

% The output, PHN the new path history, is the update of 
*PH.The decoded msg word is in the last column of PHN. (They 
sshould all have "merged " to agree on a common msg) . 

% The past histories, which are in PH, are undated on the 
sbasis of the "minimum weight". You can change this to the 
"maximum weight" if desired as indicated in the comments in 
the code if approriate for the metric of your receiver. 

6 Prior to using the function for the first symbol, you must 
initialize it with a PH matrix of dimension K by 3*Np. (I 
‘Susually use an all zero matrix). For the second and 
following symbols, PH is made the PHN from the previous 

% execution of the viterbi function. 


oO 


P=PH (22 

wt=T(:,3:3:3*2*k)’; ux=T(:,1:3:3*2*k-2)’; 

aa=D(wt(:)); bb=P(ux(:)); 

wt(:)=aa+bb; *This contains all weights(columns) for each 


state (row) 
[ a,b]=min(wt); % Use max(wt) here for maximum 
A\ayper—a 
XK (:, 1) sauce: , Soho 2) 
X(:,3)=diag(T(:,3.*b-1))}; sStores path msqword. Chg 335 
to 3.*b to keep path codewords instead of msg words. 
% Now append old paths to new paths to get survivors. 
PHN=[X PH(X(:,1),1:3*Np-3)]; 


M. MTRXV.M (OFFSET FUNCTION) 


function [a,b] =mtrxv(N,A,B) 
‘This changes matrices A and B to vectors a and b with 
soffset N. 


A=”: B=5 
a=A(:); b=B(:); 
a=a’; b=b’; 


1£ length (a) == 
a=a(li:length(a 
end 


length (b) 
) -N) ;b=b(N+1:length (b) ); 
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N. CHECK.M (ERROR CHECKER FUNCTION) 


punemton [Nb,1l,e)] =check (x,y) 
SYMBOL ERROR CHECK 
Paul H. Moose 
Naval Postgraduate School 
OS sh = S7. 


do not agree. It returns a one in e if they do not agree 
nda zero if they do agree. e 1s the error vector if x and 


6 
6 
o 
[°) 
o 
[o] 
* This m-file locates the positions in vectors x and y that 
6 
6a 
$y are binary. I is a vector of error location numbers. Nb 
ie) 
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s the sum of the elements of e. 
e=(X~=y) ; 
I=find(e) ; 
Nb=sum(e) ; 
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APPENDIX C. MSK CODES 


A.MSK.M (MSK MAIN PROGRAM) 


clear 
el 
f=100>; 
T=0.015; sSymbol Weriod 
t=0.0005;%Sample period 
ala / Modulation index 
m=5000; %tNumber of samples 
msg=random(1,m);%Transmitted binary stream 
Save msgl msg; 
msg=mapper (msg) ; 
mod output=modulato(msg,T,t,h) ; 
Save: MeGdoucl Medeouleout 
cleat Moqdmourpur 
viterbi_ path_matrix=[{1 0121 4;1 12 2 0 3]; 
coeff=(2.4841 2.2139 ]; SAWGN standart deviations 
EOnVeNE = lenge ieeer le: 
load modoutl1; 
[aa, bb] =size(mod_output) ; 
[awgn_noise] =awgn_ ch(aa,coeff(cnt) ,bb) ; 
noise=awgn_ch(aa,coeff(cnt) ,bb); 
mod_ output =mod_output+noise; 
clear noise; 
mod output=conj (mod_output’ ) ; 
Mod CUEBUE-MOGsOUEOUE (is )e: 
buf=ones (2*bb,1) ; 
mod OULDUE= | bDUE; medzouBoult |, 
clear buf; 
[bl aly =ehebyi Gadd 520. 
mod output=filter(b1,al,mod_ output) ; 
mod output-limiter (modzourpue « 
(D2yaai2'] =chebyi(4, .O25 2208); 
mod output=filter(b2,a2,mod_output) ; 
mod output=mod_ output (2*bb+1:length(mod_output),:); 
mod_output= {mod_output (8:length (mod_output) ) ;mod_output (1:7) ] ; 
mod output=reshape (mod_output,bb,aa) ; 
mod output=conj (mod output’); 
[demodulator output] =demodula(f,T,h,mod_output,t) ; 
clear mod output; 
TT=zeros (2,60); 
EOG sx — em 
D=demodulator output (x,:); 
[TT] =softv(1,1,20, FL, viterb1 jean iia masase) a, 
received signal (x) =TT(1,60); 
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end 
Clear demodulator output; 
load msgl1; 
[received_signal,msg] =mtr(19,msg, received signal) ; 
k=recelived signal~=msg; 
nofer(cnt)=sum(k) ; 
end 
Save result nofer 


B. RANDOM.M (SYMBOL SEQUENCE GENERATOR) 


function [random matrix] =random(m,n) ; 

*This M file generates m by n matrix of random binary 
values {1,0} 

mama (’ uniform’ ) ; 

rand(’seed’ ,0O); 

a=rand(m,n) ; 

b=ones (m,n) *0.5; 

random _matrix=floor (a+b) ; 


C.MAPPER.M (MAPPING FUNCTION) 


function [mapper output] =mapper (in) ; 
*Transfer 0 values to -1l 

k=in==0; 

K=-k; 

mapper output=k+in; 


D.MODULATO.M (MODULATOR FUNCTION) 
meet ion y=modulato(in, T,t,h) 


This is the implementation of the CPFSK for h=1/2. 
tetaQ=0; 


cnt=1; 
for s=l:length(in) For the number of samples 
1B aL)6\ (Veale) 
Lf Cnt== 
teta=0;%Initial theta 
else 


tetaO0O=teta0+in(cnt-1);sAccumulation of theta0d 
teta=pith*tetad; 
end 
time=0:t:T; 
time=time+(cnt-1) *T; 
th=pi*iI*h* (time/T- (cnt-1))+teta; 
mod output (s,:)=(exp(i*th)); 
cnt=cnt+1; 
end 
y=mod_output; 
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E. DEMODULA.M (DEMODULATOR FUNCTION) 


function [demed_output]=demodulartcri(t 77’, hy med .curtoue wae 
6This M file accepts modlator output and matches it with 
recreated phase samples as an input the euclidian distance 
to the euclidian distance finder and returns metric for the 
soft viterbi decoder. 
match out=[]; 

map=(-1 1]; 

(m,a] =size(mod_output) ; 
time=0:t:T; 

for sample=1:m 


mod output (sample, :)=mod_output (sample, :) *exp(j* (sample-1)*p 
Lie 
BOLT NM@any—l: 2 
xx=exp (i*time*pi*h*map(m_ary) /T) ; 
matchl=xx*cony (mod output (sample,:)’); 
match _out=[match_out matchl]; 
end 
demod_output (sample, :)=real(match_out) ; 
match out=[ ]; 
end 
R= (j*demod_output(:,1)+demod_output(:,2)); 
[demod_output] =eucdis (2,R); 


F. SOFTV.M (SOFT VITERBI DECODER) 


function PHN = softv(k,K,Np,PH,T,D) 
Soft Viterbi Decoder 
Paul H. Moose 
Univ. degli Studi di Padova 
17-05-91 


NP ol of? of? of? of? OP 


This M-file decodes k bit msgwords from 2n real 
smetrics (These may, for example, represent the "distance" 
sof the received modulation value from each of 2%n 
smodulation values.) The state transition information for a 
$2°K state trellis is in the 2°K by 3*2*k matrix T. Each of 
‘the 2°k entering paths to each state has its source state 
%(one of 2*K), path msgword (one of 2°k) and path codeword 
(one of 2°n) listed in the state row. 

The path histories are kept in matrix PH that is 2°K by 
3*Np. % The path history for each state contains source 
State, path weight and path codeword for Np previous 
states. 

The output PHN is the update of PH, the new path 
history. The decoded code word is in the last column of 


°o 


o\o o\? o\? o\? 


oP o\P 


66 


SPHN. (They should "merge". 

% The past histories are undated on the basis of the 
ltmmum Metric”. You can change this to the "maximum 
‘smetric" if desired aS indicated in the comments in the 


=D Tig 2) jeremy oie, 2) > spath weight 
Pie) =Ti7,2) gopath source state 
Pages) — (2) obaeeecode word 1(7;,3).Chg to T(j,2) for 
smsqword 
for 1=2:27k 
View (eS Seer PH (T (4 73* 1-2), 2); 
if wt < X(j,2) The < selects min metric 
Malta) =woe 
PA = Gg © 1 2a. 
PGs =a oe lose - SChg to nt, 3*1) for 
codeword. 
else 
end 


end 
* We now need to append old paths to new paths to get 
Survivors. 
PHN(},:)= (X(j,:) PH(X(j,1),1:3*Np-3)]; 
end 
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